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Two effects of the light-matter interaction are 
widely used for micromanipulation of particles 
with laser beams radiation pressure [2], orig- 
inating from a direct transfer of momentum from 
photons, and electric dipole (gradient) force act- 
ing on polarized particles [3]. In gaseous media, 
however, the heating of absorbing particles by 
light leads to much stronger radiometric forces [4] 
considered until now only as an obstacle which 
prevents a stable trapping [5]. In contrast, here 
we show that photophoretic force [6, 7] can be 
specifically tailored to trap and manipulate ab- 
sorbing particles in open air. We demonstrate ex- 
perimentally the optical guiding of clusters of car- 
bon nanoparticles [8j, [9] in the form of nanofoam 
fragments of arbitrary shape, with the size in the 
range 0.1-10 micrometers, and for laser powers 
lower than one milli Watt. The optical trap is 
created by two counter-propagating "doughnut" 
vortex beams [10, 11], and it allows simultaneous 
trapping of several particles as well as their sta- 
ble positioning and controlled guiding along the 
optical axis. Only a small fraction of operating 
power is actually absorbed because the particles 
are trapped in the region of vanishing intensity at 
the vortex core. Thus the alteration of physical 
and chemical properties of airborne particles is 
minimal in the trap, this feature is important for 
experiments with aerosols [l 2t ] . Furthermore, the 
non-contact and remote optical trapping of ab- 
sorbing aerosol particles can be applied to detect, 
monitor, and possibly reduce the exposure to en- 
gineered nanomaterials in air [13], a necessary 
tool [14] for the studies of impact of nanotech- 
nology on health [15] and environment [16]. The 
photophoretic trapping can be also employed to 
simulate, on laboratory scales, the processes stud- 
ied in atmospheric (171 . Il8| and planetary (l9l . [20( ] 
sciences. 

When a photon is absorbed by a small particle its mo- 
mentum contributes towards radiation pressure while its 
energy dissipates in heat. The latter leads to thermal 
forces [3] deliberately avoided in optical tweezers [l|: an 
optical trap utilizing a gradient force produced by a sin- 
gle strongly focused laser beam [3|. The applications of 
optical tweezers range from trapping of colloidal parti- 
cles [2l|, [2^[ and living cells [23| to manipulation of single 
molecules [24[ and atoms [25| ■ Yet these applications are 



limited by the condition of vanishing radiometric forces 
and exclude, for example, trapping of strongly absorbing 
particles in air [2l|,[26|. Indeed, if the surface of an aerosol 
particle is nonuniformly heated by an incident light, the 
gas molecules rebound off the surface with different ve- 
locities creating an integrated force on the particle, this 
effect was discovered by Ehrenhaft [|| and termed pho- 
tophoresis (PP) [1, 0]- A rough comparison [27j of the 
radiation pressure force, F a — P/c, exerted by a beam 
with power P, and the PP force, F pp = P/3w, for parti- 
cles with zero thermal conductivity [13], shows that for 
air at room temperature the later dominates by several 
orders of magnitude, F pp /F a — c/3v ~ 6 x fO 5 , here c 
is the speed of light and v is the gas molecular velocity. 
More elaborate calculations [28[ , which we provide below 
in Methods and Supplementary Notes, show that the PP 
force acting on carbon nanoclusters in air is four orders 
of magnitude larger than F a . The radiometric trapping 
of particles against gravitation, or optical levitation [27j . 
was demonstrated [29j several years prior to the first ex- 
periments on trapping with radiation pressure [2j, but 
stable all-optical trapping of particles utilizing strong PP 
forces was not realized yet. At the same time, opaque 
particles are of main concern for environmental protec- 
tion, e.g., the familiar air pollutants, such as car exhaust, 
as well as novel fabricated nanomaterials, such as carbon 
nanotubes and metallic nanoparticles. 

The major difficulty in utilizing PP force for parti- 
cle trapping is that a finite time of thermal relaxation 
within the particle is sufficient for its stochastic mo- 
tion along complex trajectories ml. Nevertheless, it has 
been experimentally observed [U, H3| that using a vor- 
tex beam with the ring-shaped transverse intensity pro- 
file [13, [HI leads to a strong two-dimensional confine- 
ment of absorbing particles on the "dark" optical axis. 
The key step forward that we make here for realizing a 
fully three-dimensional trapping is the implementation of 
the dual-beam scheme 0] but with co-rotating counter- 
propagating vortex beams [3(|, see Fig.[5]Ja). The longi- 
tudinal on-axis confinement is achieved by a balance of 
PP forces induced by two beams on the opposite sides of 
a particle while the transverse confinement by the bright 
intensity ring compensates for gravity in the horizontal 
scheme. The details of the experimental setup can be 
found in Methods Summary and Supplementary Meth- 
ods. 

To realize the PP trapping we use clusters of agglom- 
erated carbon nanoparticles [1, [9[ produced with high- 
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FIG. 1: Photophoretic trap, a, Schematic of an optical 
trap with counter-propagating and co-rotating vortex beams 
shown by surfaces at their tube-like intensity maxima. The 
focal (gray) planes of the froward (blue) and backward (red) 
beams are separated by the distance 8, for equal powers of 
two beams the trapping position is in the middle between 
two planes. Particle (green sphere) is subject to illumination 
from both sides, the geometry of the laser power flow (arrows) 
is shown with the stream-tubes, the varying width of tubes is 
proportional to the modulus of the Poynting vector, b, d, The 
shade cast by a trapped particle as seen on the optical axis (b) 
on white-light background and (d) with superimposed vortex 
beam, c, The side view of the setup with a particle trapped 
in air. A halo of the scattered light makes particle visible to 
a naked eye. 

repetition-rate laser ablation technique as described in 
Methods. The nanoclusters scatter sufficient amount of 
light to be visible by a naked eye, see Fig. [UJc). We also 
use an additional white-light source to monitor the trans- 
verse dynamics of trapped particles, as seen in Figs.[5]Jb, 
d) and corresponding Movies 1 and 2. The PP trapping is 
sufficiently robust to stabilize particle in open air for sev- 
eral minutes, as in Fig.[5{c), and we use a glass cell with 
an open top to neutralize air draughts; once captured 
in the cell the particle remains trapped for many hours 
even when the operating power is reduced below one milli 
Watt. Additional challenge of trapping aerosols [26|, in 
contrast to colloids [l|, [2l[ , is that the trapping of freely 
moving particles is passive rather than active, i.e., we 
tap the cell containing powder of fragmented nanofoam 
so that small amount is released in air around the PP 
trap, and we wait till one or several particles remain cap- 
tured stationary. In the transient regime, we observe 
fascinating scattering of many particles drawn towards 
the trapping region and competing for a stable position. 
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FIG. 2: Electron micrographs of carbon nanoclusters 
produced by laser ablation and collected from the 
optical trap, a, TEM micrograph of single nanoparticles 
collected on a holey carbon TEM grid in the laser ablation 
chamber at about 10 mm from the ablated graphite target. 
The inset shows the nanoparticle size distribution with the 
maximum at 6 nm. b, SEM image of the nanoparticle aggre- 
gates deposited on a silicon substrate approximately 200 mm 
from the laser ablation area, c, d, SEM images of samples of 
carbon nanoclusters, composite of aggregated nanoparticles, 
collected from the photophoretic trap. 



Thus we have no control over the size or the shape of 
trapped particles. 

The average size of individual nanoparticle range from 
4 nm to 8 nm as seen in the inset in Fig. [U^a). Fig- 
ure [5] also shows the Transmission Electron Microscope 
(TEM) and Scanning Electron Microscope (SEM) images 
of nanoparticles (a) near and (b) far away from the ab- 
lated graphite target, the later nanofoam is used in our 
experiments as a source of aerosol particles. The SEM 
images of nanoclusters collected directly from the trap 
are presented in Figs. [SJc, d). Among many collected 
particles of different shapes and sizes we choose to show 
one of the smallest in (c) , with the linear size of the or- 
der of 100 nm, and one of the largest in (d), spanning 
for over 10 /im. Particles smaller than 100 nm remain 
invisible in our scheme, whether trapped or not, and the 
lower limit for the size of trapped particles is unknown. 
A simple theoretical model of the PP trapping which we 
suggest and describe in Supplementary Notes II shows 
that the trapping efficiency decays rapidly for particles 
larger than the vortex ring diameter w; in our experi- 
ment this value was of the order of 10 /im which agrees 
well with the size of the largest trapped particle. 

For optical guiding we modify the dual-beam trap Q 
and include an element of optical control over the axial 
position Z of the particle. It consists of a half-wave plate 
and polarizing beam-splitter cube with low extinction ra- 
tio, approximately 1:13, see Supplementary Methods for 
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the details. Varying the tilt angle 9 of the half-wave plate 
we are able to change the power ratio e{9) = Pf/Pb 01 
the forward (P/) and backward (P&) vortex components. 
Imbalance of the powers illuminating particle from both 
sides shifts the trapping position towards a weaker beam, 
this shift is limited by the cube extinction ratio. In Sup- 
plementary Notes we provide a theoretical model which 
predicts the location Z of the trap, in a good agreement 
with experiment, as a function of the tilt angle 9, the 
size of the spherical particle a, vortex ring radius w, and 
the distance separating focal planes of two beams 5. In 
particular, we show that for a < w/2 the stationary po- 
sition Z does not depend on the particle radius a, and 
it is determined only by the separation 6 and the power 
ratio e(0). 

Experimentally we record the images of the radiation 
scattered by trapped particles in the direction perpen- 
dicular to the optical axis, the results are presented in 
Fig. [8] and Movie 3. First, we perform a static guidance 
when the particle is trapped stationary for each tilt 9 
changed with the step of 2°. The focal planes are sepa- 
rated by S = 2.0 ± 0.2 /im (s.d.e.), and we could pinpoint 
stably the particle on the beam axis anywhere within the 
distance of about 1 mm, see Fig. [5fa). It is seen that 
the position is not a strictly periodic function of 9, i.e., 
the particle does not return to the same position for a 
given value of 9. For a better comparison, we perform 
similar experiments but with the maximum magnifica- 
tion of our imaging system, the results are presented in 
Figs. [8jb, c). The half-wave plate is turned about one 
of the extrema indicated in Fig. by a dashed rect- 
angle. The uncertainty in the trapping position versus 
9 reaches values ~ 200 /im, greatly exceeding the un- 
certainty for the stationary particle, the vertical bars in 
Figs.[5Ja, b) indicate the recorded spot size ~ 10-^20 fim. 
We rule out instability of trapping because we never ob- 
served a trapped particle to change its position sponta- 
neously. The most plausible explanation for the observed 
hysteresis-type behavior is the presence of multiple traps 
produced by rather large distortion in the shape of vortex 
beams as well as misalignment of their axes. Neverthe- 
less, it is always possible to position the particle with the 
precision limited only by the resolution of our imaging 
system by appropriately adjusting e{9). 

To test the dynamic stability of the PP trapping we 
performed hand-on experiments of continuous movement 
of particles recording on-the-fly tracks with exposition 
40 ms. The half-wave plate was rotated changing e(9) 
periodically so that particle(s) were bouncing back and 
force over the distance exceeding 2 mm many times and 
reaching velocities up to 1 cm/s, here the inter-focal dis- 
tance was 8 — 7.4 ± 0.4 mm (s.d.e.). The results are 
presented in Movie 4 and Figs. [5fd-f). The position Z 
versus real time in (d) could be determined only with 
large uncertainty as the length of the tracks (vertical 
bars) recorded with finite exposition time. In this partic- 
ular experiment there were two particles trapped simul- 
taneously but they visibly separated only when moving 
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FIG. 3: Guiding of particles in air. a-c, Static guiding of 
particles: the position Z of a trapped particle measured as a 
function of the polarizer angle 9 in (a) and with higher mag- 
nification in (b). Vertical bars measure the spot size of the 
recorded particle images such as those superimposed in (c); 
corresponding data points are marked in (b) with arrows, d- 
f, Dynamic guiding of particles: the positions of two particles 
simultaneously bouncing between two extremum points of the 
trap versus time are shown in (d); black bars correspond to 
the on-the-fly tracks overlapping for both particles, such as 
the nine consecutive tracks superimposed in e. Blue and red 
bars in (d) measure the two tracks of separated particles, such 
as in the snapshot at the time t = 23 s in (f). Arrows in (e) 
and (f) indicate directions of particle propagation. 



in one direction, from the right to the left, as shown in 
the snapshot in (f). At any moment of this experiment 
the rotation could be stopped and the particles stop and 
became indistinguishable. Similarly impossible is to dis- 
tinguish two particles when their tracks overlap on the 
flight from the left to the right in (e) . Asymmetry in our 
scheme, 0.093 < e < 15.509, is attributed to loses in one 
of the interferometer arms, and the total working power 
differs by ~ 15% in two limits. Corresponding deviation 
of the Malus' law e{9) from periodic curve is clearly visi- 
ble in the map Z(9) in Figs. [5][a,b), which may facilitate 
the difference between two particles since one is blocking 
illumination of the other. 

The observations above suggest that the PP trap can 
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FIG. 4: Multiple PP trap with tilted beams, a, Volume 
plot of the longitudinal cut through the total intensity calcu- 
lated for two counter-propagating Laguerre-Gaussian beams 
LG12 tilted in the vertical direction by 0.02 rad. The yellow 
surfaces cut out the regions of small intensity where particles 
can be trapped, b, The side view and c, the front view of 
several particles simultaneously trapped with tilted beams. 



be employed for simultaneous trapping of many par- 
ticles. To elucidate this possibility we consider more 
complex multi-ring vortices created by the Laguerre- 
Gaussian beams LGi n , where the integer index n indi- 
cates, in addition to the topological charge I, the number 
of radial nodes (dark rings) in the transverse intensity 
distribution. Calculating the total intensity of the su- 
perposition of such co-rotating counter-propagating vor- 
tices we also assume a simple misalignment of their op- 
tical axes, namely a small relative tilt. This results in 
a complex light pattern with multiple minima shown in 
Fig. [7ja). Experimental results are presented in Movies 
5 and 6 with representative frames shown in Fig.[7£b,c). 
In contrast to a single trap, in a multiple trap the par- 
ticles strongly interact, as seen in Movie 5. Therefore, 
it remains unexplored whether the regular structures of 
many particles can be trapped in air. We estimate that 
the number of particles trapped in Fig. 0(a) is around 
one hundred. Combining several such beams or employ- 
ing holographic technique [![, it is possible to create a 
'web' of vortex traps in a significant volume trapping a 
large number of particles. 
SUMMARY. 

We have demonstrated a novel all-optical method for 
trapping and manipulating small absorbing particles 
in gaseous media based on photophoretic force. We 
have realized experimentally, in open air, the robust 
three-dimensional guiding of agglomerates of carbon 
nanoparticles with the size spanned for two orders of 
magnitude, from 100 nm to 10 /mi, over the distances 
of several millimeters, as well as their acceleration up 



to velocities of 1 cm/s and simultaneous trapping of 
a large number of particles. Our projections show 
that up-scaling of optical beam size will allow larger 
particles to be trapped and transported over longer 
distances, keeping trapping powers as low as few milli 
Watts. The ability of capture and controlled transport 
of nanoparticles in air by optical vortices may find 
wide applications in developing ecologically clean and 
health-safe environment for nanotechnology. The volume 
localization of nanoparticles will also allow engineering 
the appropriate control equipment for the manufactur- 
ing and micro-assembly processes. The nanoparticles 
agglomerated and collected in the non-contact and 
remote trap can be further investigated in terms of their 
chemical activity and unique toxicity (as compared to 
a bulk), important for health risk assessments. The 
outcomes are of fundamental importance for a wide 
range of other fields of science, such as interstellar dusty 
plasmas and atmospheric physics. 

METHODS SUMMARY 

Photophoretic trap. For the dual-vortex PP trap 
shown schematically in Fig. E{a), a particular care 
should be taken on the relative rotation, or topological 
charge, of two optical vortices. For paraxial beams the 
integer topological charge I determines the order of a 
phase dislocation of the complex electric field [l(| EH, 
E ~ exp(iltp + ik z z), where k z is the wave-number, 
r and tp are, respectively, the polar radius and angle 
in the transverse plane. Each time when the vortex 
beam is reflected, k z — ► — k z , the sense of rotation 
remains unchanged, I — > I, so that vortex effectively 
reverses its topological charge [3(J, defined with respect 
to the propagation direction z. Therefore, in the 
schemes based on the Fabry-Perot interferometers with 
even number of reflectors, the counter-propagating 
vortex beams, defined as E± ~ exp(iliip + ik z z) and 
E% ~ exp(— il<2Sp. — ik z z), will have opposite direction 
of rotation with the topological charges l\ = 1% = I, 
and the interference intensity pattern \E\ + -E2I 2 will be 
modulated azimuthally, i.e. ~ cos(2Zy + 2k z ). The unde- 
sired intensity minima in the transverse intensity of the 
composite trapping beam will allow particles escaping, 
and it should be avoided. Therefore, we employ a scheme 
where a single vortex beam is reflected odd number 
of times before allowing to counter-propagate itself, 
similar to a shearing interferometer. The directions of 
rotation of the initial and reflected beam coincide in 
this case, I2 = ~h, and azimuthal dependance exp(ilitp) 
is factorized in the expression of the total field, thus 
the transverse intensity distribution remains radially 
symmetric. It is also noteworthy that such a constructive 
interference effectively doubles the optical angular mo- 
mentum [30| , see the twisted power flow lines in Fig.[5](a). 

METHODS 

Synthesis of carbon nanoclusters. To produce 
samples for this study, the graphite targets were ab- 
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lated in a vacuum chamber pumped to a base pressure 
of 1(T 3 Torr and then filled with high-purity (99.995%) 
argon gas. Carbon nanofoam, consisting of agglom- 
erated nanoparticles of 4-6 nm in diameter, was syn- 
thesized using 40 W frequency doubled Nd:YVC>4 laser 
operating at wavelength A = 532 nm with repetition 
rate 1.5 MHz, pulse duration 12 ps [3l|, focused down 
to a spot size ~ 15 fim, to generate laser intensity of 
~ 10 12 W cm -2 with corresponding fluence of up to 
20 J cm~ 2 . The density of individual nanoparticles de- 
termined by the electron energy loss technique was in 
the range (1.65 1.90) g cm" 3 H [H [H[ , while the bulk 
density of the carbon nanofoam was varied in the range 
(2 20) mg cm -3 , depending on the argon pressure. 

Thermal and optical properties of nanoclusters 
(carbon nanofoam). Individual carbon nanoparticles 
have 80-90% of sp 2 bonds [H, H| and their "bulk" 
thermal and dielectric properties are similar to those 
of the graphite, namely the thermal conductivity is 
kg = 6.3 W m _1 K^ 1 (37[ and the refractive index 
N g = 1.95 + i0.66 (e g = N 2 = 3.4 + i2.6). 

We assume that nanofoam consists of a mixture of 
graphite nanospheres and air with a volume filling frac- 
tion r\ = (pf — p a )/{Pc — Pa)- Taking the average den- 
sity of the nanofoam pf — 10 mg cm -3 (author?) 
Q, of air p a — 1.29 mg cm -3 , and of a nanoparticle 
p g = 1.8 g cm -3 , we obtain 77 = 4.84 x 10~ 3 . Due to 
very low density we find, using linearized Maxwell for- 
mula kf ~ k a (l + 377), that the thermal conductivity of 
the nanofoam, kf = 0.0266 W m _1 K _1 , is mainly de- 
termined by that of air, k a = 0.0262 W m" 1 K" 1 at 
300 K (author?) [H[. 

Optical transmission measurements of the carbon 
nanofoam films with thicknesses varying in the range 
70 -j- 120 /im yield the absorption length in the nanofoam 
/ / = 35±5 pxn (s.d.e.). Corresponding imaginary part Kf 
of the nanofoam refractive index, Nf — rif + iKf, can be 
calculated as n f = X/Anlf = (1.25 ±0.15) x 10" 3 (s.d.e.). 
On the other hand, the dielectric function of low den- 
sity absorbing mixture can be found using the effective 
medium approximation, namely the linearized Maxwell 
Garnet formula [36[: £/ = e a + 3rie a (e g — e a )/(e g + 2e a ). 
Taking the dielectric constant of air s a = 1, we ob- 
tain for nanofoam Ef — 1.0079 + zO.0032 and Nf — 
1.0040 + £0.0016. The value of imaginary part of refrac- 
tive index is in good agreement with that derived from 
experimental data. 

The reflectivity of the surface of the nanofoam films, 
and thus of the nanoclusters, is very low. The measured 
diffused radiation shows no angular dependence char- 
acteristic to Lambertian source; the fraction of the in- 
cident radiation power scattered into 4tt steradian was 
(2.84 ± 0.14) x 10" 2 (s.d.e.). The power fraction of 
the reflected light at normal incidence was as low as 
(1.8 ± 0.3) x 10~ 6 (s.d.e.), similar to the value reported 
recently for the arrays of carbon nanotubes [39j] . 

Evaluation of optical forces. Theoretical treatment 
of the PP force was developed for spherical particle il- 



luminated by a plane wave. Therefore, we consider a 
sphere with typical radius a = 1 pxt\ and with thermal 
and optical parameters of the nanofoam, illuminated by 
a plane wave with characteristic intensity Jo = P/4irw 2 = 
1.1 kW cm -2 , here the vortex ring radius wq — 8.4/j,m 
and the typical power P = 0.01 W. Corresponding radi- 
ation pressure force, F a = P a /c — 4.5 x 10~ 15 N, is one 
order of magnitude larger than the gravitational force, 
F g = 4.1 x 10 -16 N, here the power absorbed by the 
particle P a = Pa s /3w%lf = 1.35 x 10" 7 W. The ab- 
sorbtion efficiency of the particle (ratio of the absorbed 
to incoming power) is of the order of a few percent, 
£ = Aa/?>l f = 0.038. 

Evaluation of the PP force is essentially more difficult. 
It involves solving consistently the electrodynamic and 
gas-kinetic problems taking into account thermophysi- 
cal, optical, and accommodation properties of an aerosol 
particle 0, [H| . We apply the limit of total accommoda- 
tion and low Knudsen number, Kn = l/a= 0.065, here 
I = 65 nm is the mean-free path of air molecules. In this 
regime the PP force is a result of the Maxwellian "ther- 
mal creep" of the gas molecules along the temperature 
gradient on the particle surface and the gas is modeled 
as continuous fluid media with boundary slip-flow condi- 
tions. The expression for the PP force is given by (see, 
e.g., formula (34) in ref. [H) 

= 9tt fi 2 a al 

pp Jl 2 Pa T(k f + 2k a y 

here the air viscosity fx a = 1.73x 10 -5 N s m -2 at temper- 
ature T — 298 K. The most important parameter which 
defines the sign and magnitude of the PP force is the 
asymmetry factor, —0.5 < J\ < 0.5, characterizing the 
angular nonuniformity of the heat sources at the particle 
surface. Depending on whether the front or back surface 
of the particle is hotter the particle will move away (pos- 
itive PP and J\ < 0) or toward the light source (negative 
PP and Ji > 0). Negative PP was observed for semi- 
transparent particles focusing illumination on the back 
side while for black body particles J\ = —0.5 and PP 
is positive. Although the nanofoam skin-depth is larger 
than typical nanocluster dimension, the thermal conduc- 
tivity is very low, thus we assume that the temperature 
changes on the particle surface are determined by the il- 
luminating intensity, similar to the black body particles. 
The factor Ji, however, takes into account corresponding 
absorption efficiency, similar to £ above, thus we adopt 
the following value: J\ = — £/2 = — 2a/3Z/. In these as- 
sumptions the PP force dominates radiation pressure by 
four orders of magnitude, F pp = 3 x 10 -11 N. 

However, the forces exerted by a light with spatially 
varying intensity, such as that of a vortex beam, can 
be essentially different. We estimate in Supplemen- 
tary Notes I that the actual radiation pressure force F a 
in our experiments is two orders of magnitude lower, 
F a = 10~ 16 N, because the particle is located on the 
dark vortex axis and absorbs much less power, P a — 
P8a 5 /15w 4 l f = 3.1 x 10~ 8 W. Thus the radiation pres- 
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sure of a "hollow" vortex beam alone might be insuffi- 
cient to counterbalance gravity. Similar reduction by a 
factor 2a 2 /w 2 is predicted for the PP force which became 
F z = 8.4 x 10~ 13 N, at least three orders of magnitude 
larger than gravitation. 

Finally, the radial trapping of particles by the trans- 
verse PP force Pr, exerted by the vortex ring, is possible 
because there is no reduction of the force magnitude, it 
scales as Pr/P z = — 8P/3a, here R is the radial displace- 
ment of the particle from on-axis stationary position, see 
Supplementary Notes II. In contrast, the transverse trap- 
ping by the radiation pressure (gradient force) requires 
strong focusing of the vortex beam by a microscope ob- 
jective with high numerical aperture, NA> 1, and the 
radial force is one order of magnitude lower than the ax- 
ial (longitudinal) force [2l| . 



Supplementary information 



Supplementary Methods: Experimental setup 

Our experimental setup is drawn in Fig. 5. The lin- 
early polarized Gaussian beam from a cw laser source 
with variable power (Verdi V5, Coherent Inc., wavelength 
A = 532nm) passes through the diffraction "fork" -type 
hologram DH where it is transformed into a Laguerre- 
Gauss vortex beam with topological charge I = 1 and 
transverse intensity pattern shown in the inset (5b) . The 
beam diameter can be varied from 2.5mm to 8mm by a 
collimator based on two lenses, LI and L2. The half- 
wave plate WP1 adjusts polarization of the input vor- 
tex beam so that it passes through the polarizing beam 
splitter BS1; the later serves as an injector of white light 
from the source WL to monitor the transverse dynamics 
of particles trapped in the glass cell C. The white light 
provides the background illumination for the images of 
the particles at the CCD2 camera after passing the notch 
filter NF which cuts off the laser radiation. 

The interferometer consists of tree reflectors - a notch 
filter NF reflecting the laser beam and passing the back- 
ground light through, a mirror M, and a polarizing beam 
splitter cube BS2. The beam splitter BS2 divides the vor- 
tex beam into two beams: the /orward-propagating beam 
(blue arrows) and the fracfciwani-propagating beam (red 
arrows) . The forward beam passes through the lenses L3 
and L4, reflects from the notch filter NF and mirror M, 
and exits trough BS2. The backward beam from BS2, 
after refection from M and NF, enters the system in the 
opposite direction and goes through L4, L3, to exit the 
interferometer reflecting from BS2. The scheme is de- 
signed so that both beams have only a single round trip, 
thus preventing unwanted interference of the beams on 
the second path. The interferometer is formed with the 
odd number of reflectors so that the axial symmetry of 
the intensity distribution is preserved for any polariza- 
tion state of counter-propagating beams, see the discus- 



sion in Methods Summary. A particle trapping volume C 
is formed between the lenses L3 and L4 and the distance 
5 between their focal planes can be varied by moving 
lens L4, see the inset (5a). The imaging camera CCD1 
collects the light scattered by the particles and monitors 
the behavior of the trapped particles in the longitudinal 
cross-section. 

The half-wave plate WP2 allows gradually change the 
tilt angle 9 of the polarization of the input optical vortex 
and thus to control the ratio e = Pf jPb of the powers 
of forward Pf and backward Pf, beams after the beam 
splitter BS2 with low extinction ratio 1:13. For full char- 
acterization of this important parameter we measured the 
powers of both beams inside the interferometer for two 
orthogonally polarized states, i.e., for 9 = and 6 = 7r/4. 
Applying Malus' law we derive the expression 



1 acos 2 26> + /3sin 2 26> 

e{9) = 5 — , 

7 1 - a cos 2 26» - f3 sin 2 20 



0.093 < e < 15.623, 
(1) 

here a = 0.928 and (3 = 0.071 are the coefficients of 
transmission through BS2 for two orthogonal linear po- 
larizations and 7 = 0.825 is the transmission coeffi- 
cient through the long arm of the interferometer taking 
into account loses of the backward beam on the mir- 
ror M and notch filter NF. The total working power 
P = Pf + Pf, inside the interferometer is less than the 
power Pj n passing from the laser onto the beam-splitter 
cube BS1 because of loses 7 < 1. It can be calculated as 
P/Pin = 7 + (1 - 7) (a cos 2 26* + /3sin 2 29) and it varies 
slightly with 9, 0.837 < P/P in < 0.987. 



Supplementary Notes 

In addition to the estimation of the magnitudes of op- 
tical forces in Methods section, based on the theory de- 
veloped for the spherical particles illuminated by a plane 
wave, here we attempt to calculate forces induced by a 
ring-shaped vortex beam. In particular, we show that the 
fraction of optical vortex power absorbed by a particle is 
two orders of magnitude lower than that from a plane 
wave and, although the full treatment of photophoresis 
problem is beyond the scope of this work, we introduce 
a simple model which allows straightforward calculation 
of the trapping position in good correspondence with ex- 
perimental data. 



Supplementary Notes I: Calculation of the 
longitudinal forces 

Radiation pressure force. The estimation of the 
radiation pressure force F a = P a /c exerted on a parti- 
cle requires calculation of the absorbed power P a . For a 
better comparison with the PP theory for a plane wave 
illumination we first introduce the characteristic inten- 
sity of such hypothetical plane wave, Iq = P/4ttw 2 (see 
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FIG. 5: Experimental setup of an optical trap for absorbing particles created by the standing wave of two 
counter-propagating vortex beams, a, Dual-beam trap with movable lens L4 adjusting the separation of focal planes S. b, 
The ring-like transverse intensity distribution of a Laguerre- Gauss vortex beam. Setup elements: DH - diffraction hologram, 
L - lenses, DP - diaphragm, WP - half-wave plates, BS - polarizing beam-splitters, WL - white light source, M - mirror, C - 
glass cell, NF - notch filter. 



Fig. EJa)), where P — J I dp is the beam power of an 
optical vortex with the intensity 

I(p, z ) = P-f e xp(--£-). (2) 

Here the ring radius w(z) — woyl-\- z 2 / z 2 with the 
beam waist wq, the diffraction length z = 2ttwq/X, and p 
is a polar radius in the plane transverse to the optical axis 
z. To calculate the absorbed power P a we use the Beer- 
Lambert law for the intensity It transmitted through the 
nanofoam film of the thickness l z , It = Ii n exp(— l z /lf), 
here 7j n is the incident intensity and If — 35 /im, see 
Methods. The power P a — J (I in — I t )dp and for a spheri- 
cal particle with the radius a we have l z — 2\J a 2 — p 2 and 
P a = 2ir /; I in {l - exp (-2Va 2 -p 2 /Z/) } pdp. For a 
plane wave I m = Iq and 

*-*-^{'-&'®}. < 3 > 

where the function f(t) = 1 — (1 + t)e _t . It is important 
for the following to note that, for small argument t -C 1, 
the function f(t) ~ t 2 (l/2 — t/3). For the representative 
particle with a = lpm (2a/lf = 0.057) we approximate 
the absorbed power as P a ~ 47r/oa 3 /3//; in terms of the 
vortex power P corresponding force can be expressed as 

P a 3 

F a ~ — - „. , for a plane wave Z in = I Q . 
c 3w z lf 



Taking typical values P — 0.01 W and w > wq — 8.4 pm 
we obtain the radiation pressure force F a = 4.5 x 10 -15 N, 
one order of magnitude larger than the gravitational 
force, F g = mg = 4.1 x 10~ 16 N, here g = 9.81 m s~ 2 
is the standard gravity and the particle mass m = 
4np f a 3 /3 = 4.2 x 10" 14 g. Note that, for small parti- 
cles a <C If, the ratio F a /F g = P/Airw 2 lfPfg does not 
depend on the particle radius a. Other parameters are 
Z = 1.1 kW cm~ 2 , P a = 1.35 x 10~ 6 W, and the absorp- 
tion efficiency £ = P a /P m = ia/3lf = 0.038, here the 
incoming power Pj n = ira 2 Io = 3.54 x 10 -5 W. 

For a spatially varying intensity /i„ = I(p, z) we as- 
sume that the spherical particle with radius a is located 
on optical axis, as shown in Fig. 2 (a), at the distance 
z = Z from the beam waist at z = 0. In our experiments 
the ring radius w > wq = 8.4 pm and zq = 837.3 pm, 
thus we can neglect the variation of the vortex intensity 
on the distances Z — a < z < Z , comparable with the 
particle size a = 1 pm. < wq -C zq, taking 

W \z-Z) = w 2 (l + Ii_lT)^ (4) 

^wl (l + = - 2 (^) - 

The incoming power is given by P ln — P f (a 2 /w 2 ) and, 
for particles smaller than the vortex ring a <C w, we 
use the expansion of the function fit) after [3] to obtain 
P in ~ Pa 4 /2w 4 = 10- 6 W. It is a factor of 2a 2 /w 2 < 
0.028 smaller than the corresponding power of a "plane- 
wave" above, P m = Pa 2 /4w 2 , because of the intensity 
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FIG. 6: To the calculation of the longitudinal forces, a, 

Red-shaded is the part of a vortex beam illuminating particle 
(green sphere), b, Transfer of a momentum (red arrow) from a 
gas molecule to a particle; the illuminated side of the particle 
is a hemisphere n < 6 < n/2. c, The PP force [5] versus axial 
distance z, the bottom frame shows the ring radius w(z) of 
the diffracting vortex beam, d, Amplitude (f(a 2 /wo), top) 
and the full width at half maximum (FWHM, bottom) of the 
PP force curves shown in (c) . The red and blue curves in (d) 
show the lowest order approximations of [S] f(t) ~ t 2 /2 (red) 
as in[6] and taking into account next term f(t) ~ t 2 (l/2 — t/3) 
(blue). 



zero at the vortex origin, see the red-shaded region in 
Fig. 2 (a). Similarly, the absorbed power is two orders 
of magnitude lower, P a = 3.1 x 10~ 8 W, than that for a 
plane wave above (the absorption efficiency is of course 
the same £ = 16a/15Z/ = 0.031); to the leading order 
of the small parameter a/w -C 1 the radiation pressure 
force is given by 



F, 



I5w 4 l 



for a vortex beam 7;, 



/ 



I(p,Z). 



Therefore, the force F a = 10~ 16 N exerted by a vortex 
beam is insufficient to compensate gravity for a chosen 
level of optical power. Note, however, that the ratio 
F a /F g scales as a 2 (limited by the assumptions above 
a -C w and a -C If), so that the radiation pressure 
can overcome gravity for larger particles and for higher 
optical power. 

Photophoretic force. As we briefly discussed in 
Methods, the calculation of the PP force involves solving 
consistently two problems, the electrodynamic radiation 
problem including absorption and re-emission, and the 
gas-kinetic equations taking into account thermophysi- 
cal and accommodation properties of an aerosol particle. 
Most of the works on PP employ Mie scattering theory 
to solve the first problem, therefore, the available solu- 
tions, such as Eq. (1) in Methods, were derived for a 
spherical particle illuminated by a plane wave. Here we 
attempt to overcome this limitation and take into ac- 
count the spatially varying vortex intensity. In general, 
the solution to the electrodynamic problem with non- 



uniform illumination can be obtained, for spherical par- 
ticles with given optical properties, as a superposition 
of known Mie solutions for plane waves. However, the 
nanoclusters have random shape and they are composite 
of carbon nanoparticles, thus the optical properties are 
determined by a complex multiple scattering and can be 
only estimated using rough effective medium approach, 
see Methods. Therefore, instead of involving calculations 
with Mie solutions we adopt a somewhat intuitive theo- 
retical approach. 

As we demonstrated above the amount of optical power 
absorbed by a nanocluster is very low because of a large 
absorption depth in nanofoam, therefore, we can neglect 
the distortion of vortex beam intensity by the particle. 
Then, because of low thermal conductivity we assume 
that the temperature changes on the particle surface are 
determined by the illuminating intensity 7, similar to the 
black body particles. Finally, we introduce the linear mo- 
mentum flux, \M\ = kI, transferred from gas molecules 
to the particle, as shown in Fig. \T[h). The phenomeno- 
logical coefficient k has the dimension of inverse velocity 
and absorbs thermal and optical parameters of the gas 
and the particle. By definition, the PP force is given as 
an integral over the momentum flux density, F — J MdS, 
here dS — 2ira 2 sin 0d9 is particle surface element taking 
into account cylindrical symmetry of the problem. The 
longitudinal PP force F z is given by (see Fig. E[b)) 



M z dS, M z — — k cos 7(5+) , 



(5) 



here 5+ is illuminated hemisphere, tt/2 < 8 < tt. For il- 
lumination by a plane wave with characteristic intensity 
7o we obtain F z = K7ra 2 7o. Comparison with the expres- 
sion for PP force, F pp = 3 x 10~ n N, derived in Methods, 
allows us to evaluate 



K = — J l 



9^ 



2ap a T(k f + 2k a ) 



_ = 8 5 x iq-7 £ 

lfp a T(kj + 2k a ) m' 

here the last expression is obtained with Ji = — £/2 = 
—2a/3lf (see Methods), derived above for small parti- 
cles a < I/, the thermal conductivity of the nanofoam 
kf = 0.0266 W m _1 K _1 and the parameters of air 
are: viscosity fi a — 1.73 x 10 -5 N s m -2 , temperature 
T = 298 K, mass density p a — 1.29 mg cm -3 , and ther- 
mal conductivity k a = 0.0262 W m _1 K _1 . Note that 
the final expression for k does not depend on particle 
radius and optical power, it characterizes the transfer of 
momentum from air molecules to the material with given 
properties. 

For a vortex beam the intensity 7(5+) in O is taken 
from [5] with the approximation [5] Integration in [5] gives 
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FIG. 7: To the calculation of transverse PP force, a, 

Geometry of the problem; arrows show the transverse projec- 
tion of the momentum flux density M± transferred to the par- 
ticle (green sphere); corresponding amplitude |Afj_| = \M P \ is 
color-coded in b. The background in pictures (a, b) is the 
grey-coded intensity of a vortex with the ring radius w. c, 
d, The magnitude of the transverse force [8] versus two key 
parameters, R/w and a/w. 



with the function f(t) defined in [3] For small particles 
a <C w the expression can be greatly simplified, 



Pa 4 

b z ~ k 7, 

2 w 4 ' 



(6) 



and it differs from the solution for plane wave illumina- 
tion above, F z = kPcl 2 /Aw 2 , by a factor of 2a 2 /w 4 < 
0.028, exactly the same reduction as for radiation pres- 
sure force. Therefore, the actual force exerted by a vortex 
beam on our representative spherical nanocluster with 
a = 1 /im is F z = 8.4 x lO" 13 N. 

The PP force [5] can be measured in units of kP and it 
is visualized in Fig. [Bfc, d). The bell-shaped curves in 
(c) appear because of vortex diffraction, w(z), shown in 
the bottom panel for better comparison. As seen in (d) 
the amplitude of the force saturates for large particles 
a > 2wo so that, for the particles with a < 0.5wq, the 
approximation [5] is sufficiently accurate while for larger 
particles the next term in the expansion of f(t) can be 
taken to improve the results. The full width at half maxi- 
mum (F WHM) of the curves in (c) has weak dependance 
on the size parameter a/wo as seen in the bottom panel 
in (d). 



Supplementary Notes II: Calculation of transverse 
PP force 



To estimate the geometry and the magnitude of PP 
force in transverse plane F± and evaluate its trapping ef- 
ficiency we adopt the same assumptions as for the calcu- 
lation of the longitudinal force above, namely we assume 



a spherical strongly absorbing particle with the radius a 
much smaller than the diffraction length zq so that [5] ap- 
plies. Furthermore, from the geometry of the momentum 
transfer depicted in Fig. [f^b) we can see that the modu- 
lus of the area density of the transverse momentum flux 
is given by the projection (cf. [S]) 



K sim 



while its direction is towards the particle origin. Since 
sin# = p/a, where p is the polar radius in the transverse 
cross-section of the particle, we derive: 



Mx = -k-I(S+), 
a 



F, 



M \dS 



s+ 



I(S)dS. 



We assume that the particle is shifted by a vector R 
from the optical axis, then the distribution of the momen- 
tum flux density is mirror-symmetric with respect to R as 
shown in Fig. U^sl) . There are two zeros in the distribu- 
tion of the amplitude \M± \ as seen in (b), one zero at the 
origin of the vortex intensity ring and another at the ori- 
gin of a particle because of zero projection at 9 — n. The 
combination of two minima produce highly asymmetric 
distribution along R which is the source of the PP force 
returning particle to the vortex origin. From the mirror 
symmetry follows that the transverse force is antiparal- 
lel R and directed towards vortex axis, F± = ^Fr and 
Fr < 0. Therefore, for simplicity and without loss of 
generality we can choose, e.g., the cartesian projections 



R x — and R y = R, 



F R = — 



kP 



y x 2 + {y + Rf 



x 2 -\-y 2 <a 



(7) 



x exp 



(y + Rf 



w 2 {Z) 



adxdy 



Transverse force is calculated as a function of two key 
parameters, R/w and a/w, the results are presented in 
Fig- EI C 7 d). Several conclusions can be readily made. 
First, there is a sharp boundary between the trapping 
region inside the beam (R < w) and the repulsion of par- 
ticles far from optical axis (R > w); the unstable equi- 
librium (F p = 0) is located approximately at the ring 
radius R < w. Second, while the small particles are al- 
ways trapped in the vicinity of optical axis, for particles 
larger than the ring a > w the trapping efficiency rapidly 
vanishes; the maximal trapping force is achieved for par- 
ticles with a ~ 0.85w. Finally, the remarkable scaling 
property of [8] leads to a somewhat counterintuitive re- 
sult: the transverse force on a large particle (a > 0.85w) 
can be stronger for lower vortex intensity when particle 
moves away from the beam focus. Indeed, the radius 
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w(Z) increases because of the diffraction and the param- 
eter a/w decreases, thus the particle effectively follows 
the curves in (d) from right to left passing by the force 
maximum. The reason for this unexpected result, as seen 
in (a, b), is the interplay of geometrical contributions to 
the momentum flux versus spatial structure of the vortex 
intensity "weighting" momentum flux density in[5J 

Transverse force vanishes at the vortex axis, R — > 0, as 
should be expected from symmetry considerations. Close 
to this equilibrium position and for small particles, a + 
flCwo, we can linearize [8] and obtain 



4 

kP- 



Ra 3 
3 w 4 (Z) 



4 Ra 3 
< kP- 



— ■ (8) 



3 w. 



In our horizontal scheme the transverse PP force should 
be compared to the gravitational force, F g = 4npf a 3 g/3, 
here g = 9.81 m s~ 2 is the standard gravity and pf — 
10 mg cm~ 3 is the characteristic mass density of the 
nanofoam, see Methods. Because the PP force |S] scales 
as a 3 , similar to F g , it is possible to increase particle size 
and still balance gravitation without increasing the laser 
power P. 



Supplementary Notes III: Calculation of the trap 
position 

Based on the setup in Fig. [5] we have two principal de- 
grees of freedom for spatial manipulation of trapped par- 
ticles, apart from the total power P of the beams which 
defines the magnitude of the PP force. One is the dis- 
tance 8 between focal spots of two counter-propagating 
beams, adjusted by the position of the lens L4. Another 
is the tilt 9 of the half- wave plate WP2 defining the power 
ratio of two beams e(6) in [TJ Each beam applies a PP 
force given by [5] being a function of the position of the 
particle, here z, relative to the beam focal plane z = 0. 
For two beams with focuses separated by the distance 8 
we introduce the origin of the coordinates in the middle 
between focal planes located at z = ±6/2, see Fig. EJa). 
The total longitudinal force on the particle is then given 
by the sum, F tot = F z (z - 8/2; P f ) -F z (z + 8/s; P b ), so 
that using [5] we obtain 



Ftot (z) 




(9) 



here Pf and Pb are the powers of the forward and back- 
ward beams, respectively. Stationary points z = Z can 
be found from the condition F to t{Z) = 0. As we demon- 
strated above, for the particles with 2a < wq < w(z) 
the linearized expression for the force [HI provides good 
approximation. In this limit the stationary point (trap 
position Z) does not depend on the particle radius a and 



it is given by the solutions to a corresponding quadratic 
equation 



+ 8e+l 

Z ± = — ± 

2e- 1 



S 2 ( 



(6-1)2 



-0- 



(10) 



here e 2 (9) — e(9) from[T] However, the stationary point 
can be stable or unstable. Linearizing [TU] in the vicinity 
of a stationary point, z = Z + £ with |£| -C \Z\, we obtain 
F to t(Z + £) ~ £ (dFt t/dz) \ z -z- Therefore, the station- 
ary point Z is stable if (dF to t/dz) \ z= z < (the force 
is returning), and it is unstable otherwise. Calculation 
of the derivative of the total force [10] gives the following 
result: 



rfFto. 

dz 



5=Z± 



±8 Z ± C\ here C (8, Z±) > 0. (11) 



Therefore, there is one stable and one unstable station- 
ary point, see Fig. [UJa), and the stability is defined by 
the sign of the root in \W\ as well as the sign of the 
inter-focal distance 8. For equal powers of two beams, 
P f = P b and e = 1 at 9 = tt(0.134, 0.366, 0.634, 0.866), 
the force derivative is proportional to —8, instead of flTT 
and the trapping position Z — is stable only for <5 > 0, 
when focal planes of two beams are separated as shown 
in Fig. EJa). The case 8 < corresponds to the anti- 
separated focuses; although there are stable trapping po- 
sitions, they are located outside the inter-focal region 
where the transverse force rapidly decreases, see Sup- 
plementary Notes II above. Furthermore, the trapping is 
possible (the roots in[10]are real) only if the separation of 
focuses 8 exceeds some minimal value, \8\ > zq\e — \\/y/e. 
In the opposite case the real roots disappear and there 
is no stationary trapping on axis. Similarly, if we use 
an ideal polarizing beam splitter BS2 in our scheme, so 
that a = 1 and (3 = in [TJ the power ratio would be- 
come unbounded e = tan 29/^/ and, for any 8, the trap 
position with disappear (move to infinity Z — > oo) for 
some value of 9. It could still be used for stationary 
guiding (positioning) of particles, as in the manuscript 
figure 3(a-c), but the dynamical bouncing of particles, 
as in the manuscript figure 3(d-f), will be impossible. In 
contrast, for our "unperfect" beam-splitter BS2 we have 
bounded power ratio, 0.304 < e(9) < 3.953, and there is 
the stationary solution [TOl for any < 9 < 2ir (full turn 
of the input polarizer) if 8 > 1.485zq- 

Our objective here is to find the domain L of stable 
trapping along optical axis as a function of 8 for the given 
characteristic of a beam-splitter [TJ For separated focal 
planes <5 > this domain includes a middle point Z = 
for equal powers and it is given by L = Z~(e max ) — 
Z+( ^min); here 

(i) the inter-focal distance 8 < 1.261 zq: stable trapping 
is limited by a finite domain of 9, for a given 8 the limits 
satisfy e ma x,mi„ = 1 + A/2 ± A^l + A 2 /4, here A = 
8/zq, see an example in Fig. [S^b). The domain is L — 
^J8 2 + 4zfi so that, for 8 — > and 1 ± 6/zo, 
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FIG. 8: To the calculation of the trap position, a, The 

total longitudinal forces (black solid lines) for equal (e = 1) 
and unequal (e — 4) powers at the inter-focal distance S — 
zq. Corresponding stable (S) and unstable (U) positions of 
the trap are indicated in the Z(9) diagram in (b). Similar 
diagrams (c) and (d) show stable (red) and unstable (blue) 
roots Z(8) for different values of 5: solid lines for 8 > and 
dashed lines for 5 < 0. The calculated position of the trap (red 
lines) in (e) and (f) is compared to experimental data (black 
bars) from the manuscript figures 3(a) and 3(d), respectively. 



we formally obtain the minimal value of the domain L — > 
2zq. However, at the limit (5 = and e — 1 the force is 
zero everywhere, so that close to this limit the magnitude 
offeree is infinitcsimally small. The domain growths with 
5 — ► 1.261 zo and L m ax = 2.364 zq. 

(ii) the inter-focal distance 1.261 < S/zq < 1.485: here 
e min = 0.304 and e max = 1 + A 2 /2 + A^l + A 2 /4. 
The domain L decreases monotonically from 2.364^0 to 



2.015 Zo, the trapping is continuous on a limited interval 
of 0, see Fig.lSJc). 

(Hi) the inter-focal distance 1.485 zq < 5: trapping 
is continuous for £ [0, 27r) with e m i n = 0.304 and 
£max = 3.953, see Fig. [8jd). The trapping domain, 
L/z Q = 1.7762A - V0.4534A 2 - 1 - V0.6290A 2 - 1, 
first decreases to reach its minimum L m i n = 1.3737 zq at 
5 = 2.4351 zq, then diverges with increasing S. 

We compare theoretical calculations with experimen- 
tal results in Fig. [5]^e,f). Both experiments were done 
for continuous trapping in the whole domain 6 € [0, 2ir), 
described in the case (Hi) above with the corresponding 
roots shown in Fig. [5£d.) . For the static guiding (posi- 
tioning) in Fig. [Hie), as in the manuscript figure 3(a), 
the inter-focal distance S = 2.04 mm corresponds to the 
minimal value of the trapping domain L — 1.15 mm, 
slightly larger than the observed guiding over the dis- 
tance ~ 1 mm. However, the slopes of the curve in 
Fig. [5^e) agree well with experimental data points. It is 
not the case in Fig. (8jf) , where the results of dynamical 
guiding, from the manuscript figure 3(d), are compared 
to the calculations of stationary trap position. Indeed, 
the slopes of the experimental dataset are much steeper 
than that of the solid theoretical curve because the slope 
of the curve Z(t) corresponds to the particle velocity, 
the latter reaches values up to ~ 1 cm/s, determined by 
measuring the longest particle tracks (vertical bars in (f)) 
~ 400 /im recorded with exposition 40 ms. Furthermore, 
the extrema of the experimental dataset are flatter than 
that of the theoretical curve, this difference is a direct 
manifestation of the inertia of trapped particles moving 
along the optical axis and spending more time close to the 
points of return. The discontinuity visible in the middle 
of the dataset in (f) is because the rotation of the half- 
wave plate was not monotonic (simply by hand). More 
importantly, the theoretical predictions on the trapping 
domain L ~ 2.4 mm, calculated for the inter-focal dis- 
tance 8 ~ 7.3 mm used in experiment, agree well with 
the limits of particle motion L cxp ~ 2.2 mm. We con- 
clude that relatively simple theoretical model developed 
here can be used to estimate and predict the domain of 
the PP trapping. 
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